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We present an extensive analysis of new light and radial-velocity (RV) curves, as well as high- 
quality broadening-function (BF) profiles of twelve binary systems for which a contact configuration 
with large temperature differencies between components has been reported in the literature. We 
find that six systems (V1010 Oph, WZ Cyg, VV Cet, DO Cas, FS Lup, V747 Cen) have near- 
contact configurations. For the remaining systems (CX Vir, FT Lup, BV Eri, FO Hya, CN And, 
BX And), our solutions of the new observations once again converge in a contact configuration with 
large temperature differencies between the components. However, the bright regions discovered in 
the BFs for V747 Cen, CX Vir, FT Lup, BV Eri, FO Hya, and CN And, and further attributed to hot 
spots, shed new light on the physical processes taking place between the components and imply the 
possibility that the contact configurations obtained from light- and RV-curve modelling are a spurious 
result. 
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Theoretical investigation of W UMa-type binary stars began with the work of 
Kuiper (1941), who found that contact binary stars with components belonging to 
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the Main Sequence (MS) cannot exist in an equilibrium state (a result known as the 
Kuiper paradox). Lucy (1968) has shown that Kuiper's conclusion is valid only for 
contact binary stars with radiative common envelopes and unequal-mass compo- 
nents: contact binaries consisting of MS stars with common convective envelopes 
can exist in hydrodynamical and thermal equilibrium states if CNO and p-p cycles 
occur in the cores of both components. According to Lucy's common-envelope 
model, both stars in a contact binary system should have equal effective tempera- 
tures. Although this result resolved the Kuiper paradox, Lucy's theoretical mod- 
els of contact binary stars with components differing in their structure because of 
CNO and p-p cycle reactions covered only the blue end of the period-color diagram 
(Eggen, 1967). 

The small masses of primary components obtained from combined radial- velo- 
city (RV) and light-curve modelling (Lucy, 1973) prompted Lucy (1976) and Flan- 
nery (1976) to consider a situation with both solar-type components harbouring 
p-p reactions. It turned out that in such a situation it is not possible to attain a 
thermally stable equilibrium. Instead, in a thermal time scale of 10 milion years, 
systems will oscillate between contact and semi-detached configurations. Robert- 
son & Eggleton (1977) performed more realistic computations, taking into account 
an angular-momentum loss during a system evolution. As a result, they obtained 
models of contact-binary systems with physical parameters of components closely 
resembling those obtained from light- and RV-curve modelling. 

According to the leading work on the Thermal Relaxation Oscillation (TRO) 
theory (Lucy, 1976; Flannery 1976), the maximum temperature difference (AT) 
between components of solar-type contact-binary system should not exceed 880- 
1000 K. Similar values have been recently obtained by Li et al. (2004a, 2004b, 
2005) and Kahler (2004). The marginal-contact configurations and temperature 
differences between the components in the 300-700 K range obtained for RW PsA, 
AK Her, and W Crv by Lucy & Wilson (1979) have confirmed the theoretical find- 
ings. 

Kaluzny (1983, 1985, 1986), Hilditch et al. (1984), and Hilditch & King (1988) 
have found other contact systems exhibiting large AT between their components. 
For some of these systems, light-curve modelling resulted in temperature differ- 
ences that were considerably larger than the theoretical limit of about 1000 K. 
Rucinski (1986) and Eggleton (1996) proposed that indications of contact binaries 
with such large temperature differencies are in fact artefacts caused by the appli- 
cation of simple Roche lobe-based light curve synthetic models, with disregard for 
intensive accretion processes expected to take place between components during 
short and violent phases when a contact configuration is broken or re-established. 

In Table 1, we present a sample of close binary stars for which at least one 
light-curve solution was previously obtained using a light-curve synthesis program 
based on the Roche model and yielded a contact configuration with a temperature 
difference between the components larger than 1000 K. However, as indicated in 
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column 7 of Table 1, almost all the models are based on photometric mass ratios, 
a limitation which significantly decreases the reliability of the obtained configura- 
tions. Closer inspection reveals that in many cases an O'Connel effect (i.e., occur- 
rence of unequal maxima, most probably due to a combination of effects introduced 
by spots on one or both components) has not been taken into account during the 
light-curve modelling. Additionally, closer inspection reveals the applied albedo 
and gravity-darkening coefficients to have been far from their theoretical values. 
Moreover, the scatter visible in the light curves was often too large to permit a 
unique solution. We refer an interested reader to references given in the last col- 
umn of Table 1 for more historical information about the targets investigated in this 
paper. 

All the above has led us to the conclusion that in order to get unique con- 
figurations for these systems, one must obtain high quality light and RV curves 
(Section 2) and must perform advanced modelling on this data (Section 3). Re- 
sults obtained from the modelling of individual systems are described in Sections 4 
and 5. Physical parameters of these systems are presented in Section 6. In Sec- 
tions 4 and 5 we also present a discovery of bright regions (hot spots) on secondary 
components, made possible by the outstanding quality of the broadening functions 
(BFs) obtained for some systems. In Section 7, we discuss the implications of this 
discovery for results obtained classically (i.e., obtained from light- and RV-curve 
fits within the Roche-lobe model). 

2. Observations and data reduction 

2.1. Spectroscopy 

Medium-resolution spectra of FS Lup were obtained at the South African As- 
tronomical Observatory (SAAO) in the second half of April 2006, with the 1.9 m 
Radcliffe telescope equipped with a slit spectrograph and the grating "no.4" (1200 
1/mm). With the 1.2 arcsec-wide slit projected on 2.2 pixels of a CCD camera, this 
results in effective resolving power R = 8 000 . The star was observed in two spec- 
tral regions: during the first night the spectral window was centered at 4400A, and 
during all other nights at 5250A. 

FT Lup, CX Vir, and V747 Cen were observed in service mode at ESO-La 
Silla during six nights of April, May, and June 2006. In this facility the 2.2 m 
MPG telescope equipped with the FEROS spectrograph provided high resolution 
(R = 48000) and very good signal-to-noise spectra. 

The targets reachable from the northern hemisphere (WZ Cyg, VV Cet, BX And, 
DO Cas, and BV Eri) were observed at the David Dunlap Observatory (DDO), a 
suburban facility then owned by the University of Toronto, during November and 
December of 2006 and February of 2007, with the DDO 1.9 m telescope and a 
slit spectrograph. Most of the data were obtained with the 2160 1/mm grating, but 
because of the relatively faint magnitude of VV Cet and WZ Cyg, during nights 
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Table 1 

The sample of contact binary systems with large temperature differencies between the components. 
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The line in the middle separates systems with radiative envelopes (above) from systems with 
convective envelopes (below). The column denoted "conf." gives the configuration derived for 
the system in question: C for contact; m-c for marginal contact; SD1 and SD2 for semi-detached 
systems, where it is, respectively, the 1st or the 2nd component that fills its Roche lobe; and n-c 
for near-contact. The quantity / is the fill-out parameter, given only for contact and marginal- 
contact systems in cases in which we had the necessary data to make a computation through the 
equation /=(£! — £2^ )/(Q.i l — Q.i 2 ) (where Q,^ and Q.i 2 are the Roche potentials at the L\ 
and L2 points, respectively, Q. being a common surface potential obtained from the light-curve 
modelling). We also show the mass ratios q obtained from the light-curve solution, adding an 
asterisk to mark cases in which a spectroscopic mass ratio was used. 



of poor seeing (a few arcsec), the observations were made with a lower-resolution 
(1200 1/mm) grating. With the 1.5-arcsec slit projected on 2.5 pixels, the effective 
resolving power was R = 19000 and R = 1 1 000 for the 2160 1/mm and 1200 1/mm 
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gratings, respectively. 

FO Hya was observed at the Las Campanas Observatory (LCO) in February 
2009 with the 2.5 m Du Pont telescope, equipped with the classical slit echelle 
spectrograph. With the 1-arcsec slit, this instrument assured spectral resolution of 
R = 40000. 

The data gathered at the ESO-La Silla observatory were reduced in the usual 
way with FEROS-DRS under ESO-MIDAS. The data obtained at SAAO, DDO, 
and LCO were reduced with IRAF. 

For RV measurements, we applied the BF method (Rucinski, 1992, 1999, 2002, 
and Pribulla et al., 2009; the 2009 publication has the fullest list of references). BFs 
were calculated on the basis of RV standards from the list of Stefanik et al. (1999), 
observed in all the observatories on the same nights as the main targets, and on 
the basis of synthetic spectra from the SPECTRUM software (Gray, 2001), using 
Kurucz's (1993) atmosphere models for MS stars with solar abundances. 

RVs of light centroids were measured by means of synthetic rotational profiles 
for spherical stars. This choice assured reasonable model-independence, while still 
approximating very accurately the observed profiles of elliptically disturbed stars. 
As primary and secondary components of our targets have very different spectral 
types, the RVs of primary components were obtained from BF profiles calculated 
using the early-type RV standards, while RVs of secondary components were ob- 
tained from BF profiles calculated using late-type RV standards (Table 2). We 
found that in most cases this approach minimized scatter of RV curves and yielded 
almost identical (within 1 —2 km/s) mean system velocities v y for the primary 
and secondary components. The BF method gave larger semi-amplitudes of the RV 
variations (Kj) than did the cross-correlation function (CCF) method previously ap- 
plied for RV measurements of CX Vir (Hilditch & King, 1988), FT Lup (Hilditch et 
al., 1984) and BX And (Bell et al, 1990), thereby fully confirming the predictions 
of Rucinski (1999) that there would be a tendency for systematic reduction of Ki 
under the CCF approach. 

Table 3 contains preliminary orbital parameters obtained from the new RV 
curves by fitting the sine-curve to the data in the phase ranges 0.15-0.35 and 0.65- 
0.85. 

2.2. Photometry 

New high-quality photometric data for targets reachable from the northern hemi- 
sphere were taken during the three consecutive northern autumn-winter seasons 
2004/2005, 2005/2006, and 2006/2007. Most of the data were obtained with the 
Carl-Zeiss 50 cm telescope at the Fort Skala Astronomical Observatory of the 
Jagiellonian University. This telescope was equipped with the Photometries S300 
CCD camera with the SITe SI003B, 1024x1024 pixels chip, carrying a set of broad- 
band BVRI (Bessell) filters. With a focal length of 6.7 m, the field of view was 
12 arcmin x 12 arcmin. In this observatory, we took complete light curves for 
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Table2 

Spectral regions and RV-standard stars used for the RV determination. 



system 


used spectral range [A] 


RV standards 


FSLup 


3920-4750 & 4880-5640 


HD 102870 (F8V) & HD 92588 (KHV) 


FTLup 


4378-5743 (Hp excluded) 


F2V & K3V 


CX Vir 


4378-5743 (Hp excluded) 


F5V & K0V + K3V 


V747 Cen 


4905-5743 


A5V & K1V 


FOHya 


4910-5560 


F0V & K0V 


V1010 Oph 


5074-5305 


HD 128167 (F2V) & HD 144571(G8V) 


BVEri 


5074-5305 


HD 128167 (F2V) & HD 65588 (G8V) 


BX And 


5074-5305 


HD 128167 (F2V) & HD 3765 (K2V) 


DO Cas 


5074-5305 


A5V & K1V 


VV Cet 


5074-5305 


HD 222738 (F7V) 


WZCyg 


5074-5305 


HD 222738 (F7V) & HD 3765 (K2V) 



Table3 

Preliminary orbital parameters obtained by sine-curve fitting to spectroscopic observations in the 
phase ranges: 0.15-0.35 and 0.65-0.85. 
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DO Cas 


78.85 
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VV Cet 
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WZCyg 


130.91 
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1.05 


205.11 


± 


3.54 


0.98 


± 


1.52 
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± 


0.016 



WZ Cyg, BX And, GW Tau, and BL And. 

DO Cas, VV Cet, and CN And were observed at the Mount Suhora Observatory 
of the Cracow Pedagogical University using the Carl-Zeiss 60 cm telescope. The 
first two stars were observed with a three-channel photometer equipped with the 
Hamamatsu R2949 photomultipliers and broadband Johnson BVRI filters. CN And 
was observed using the SBIG ST10XME CCD camera (containing the KODAK 
KAF-3200E/ME chip) and wideband BVRI (Bessell) filters. 

CX Vir and FT Lup were observed in April 2005 at the SAAO with the 50 cm 
Boiler & Chivens telescope and a single-channel modular photometer equipped 
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with the Hamamatsu R943-02 GaAs photomultiplier and B V (RI) C Johnson-Cousins 
filters. In March 2006, at the same observatory, we observed V747 Cen with the 
75 cm Grubb-Parsons telescope, equipped with the UCT CCD camera with the 
EEV 576x420 pixel chip (binned 2x2) and Johnson BVRI filters. The same tele- 
scope, but equipped with the single-channel UCT Photometer with the Hamamatsu 
RCA 31034A GaAs photomultiplier and Johnson-Cousins BV(RI) C filters, was 
used during April and May 2006 for observations of FS Lup and FO Hya. 

Data gathered taken with CCD cameras were reduced in the standard way 
through ESO-MIDAS scripts. Aperture photometry with variable aperture size 
was performed through Kopacki's programs (private communication), based on 
the DAOPHOT II package (Stetson, 1987). Data obtained with single-channel 
photometers were dead-time corrected using the procedure given in the XLUCY 
manual (Balona, 2000). Sky-background and comparison star interpolations were 
applied with our own software. The data taken with the three-channel photometer at 
the Mount Suhora Observatory were reduced using dedicated computer programs 
developed by M. Drozdz (private communication). Standard corrections for differ- 
ential extinction and colour extinction were applied to the new light curves. 

Light and RV curves were phased with the most up-to-date elements deter- 
mined by Kreiner (2005-2008, private communication). These ephemerides were 
obtained using the minima times included already in the Kreiner (2004) database, 
but supplemented with several new times of minima determined from our own ob- 
servations. The comparison and check stars, as well as elements applied for phases 
calculation, are listed in Table 4. 

As attempts to take RV curves of GW Tau and BL And failed in bad weather, 
these targets were excluded from the current investigation. 



Table4 

Comparison and check stars applied for differential photometry. The two right panels contain 
elements used for phases calculation. 



system 



comparison star(s) check star(s) 



M [HJD] 



BX And 
WZCyg 
VV Cet 
CN And 
FO Hya 
CX Vir 
FT Lup 
FS Lup 
V747 Cen 



mean of 3 
mean of 3 
GSC 4674 
GSC 2787 
GSC 6049 
GSC 6137 
GSC 7828 
GSC 8292 
GSC 8261 



stars 
stars 
0662 
1803 
0074 
0576 
1478 
1678 
1863 



internal comp. 

internal comp. 

no check star 
GSC 2787 1891 

no check star 
GSC 6138 0644 
GSC 7829 0247 
GSC 8292 0914 
GSC 8261 0719 



2453350.23061 
2454019.26728 
2453675.39228 
2452500.12040 
2453855.35760 
2453470.50082 
2453478.64697 
2452500.08800 
2453811.34417 



0.61011240 
0.58446762 
0.52239410 
0.46279081 
0.46955663 
0.74608003 
0.47008002 
0.38139970 
0.53719482 
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3. Light- and RV-curve modelling: the method 



Light- and RV-curve modelling was performed through the Wilson-Devinney 
(WD) code (Wilson, 1996), with a Monte Carlo search method applied to locate a 
global minimum in parameter space (Zola et al., 2010). An iterative procedure of 
light- and RV-curve modelling, introduced by Baran et al. (2004), was used during 
our computations. Errors of adjusted parameters were calculated from the search 
matrix at confidence level of 90%, as described by Kreiner et al. (2003). The values 
of a reduced (weighted) % 2 were used to determine the quality of the obtained fits. 

The effective-temperature- versus-spectral-type calibration of Harmanec (1988) 
was used to determine the effective temperatures of the primary components ( 7^- im ) 
according to their spectral types. (Table 5 gives details.) This parameter, like 
the spectroscopic mass ratio ^ spec , was not adjusted during light curve modelling. 
The albedo A and gravity-darkening coefficients g of both stars were assumed 
at the theoretical values, i.e. A = 1, g = 1 for stars with radiative envelopes 
(r eff > 7200 K) (von Ziepel, 1924), and A = 0.5, g = 0.32 for stars with con- 
vective envelopes (T e s < 7200 K) (Lucy 1967, Rucinski 1969). The square-root 
limb darkening law and the coefficients published by Diaz-Cordoves et al. (1995) 
and Claret et al. (1995) were assumed as a function of an effective temperature of 
a star and the effective wavelength of the filter. As the choice between the detailed 
and simple reflection-effect treatment turned out to have no practical impact on the 
obtained results, especially on the obtained configurations, we fitted our data under 
the assumption of a simple reflection model (mref =1) with a single reflection 
(nref = 1). 

During light curve modelling, we adjusted the orbital inclination i, the sec- 
ondary-compontent effective temperature r e s | c , the Roche potentials £2 and the 
luminosities L of both components, the phase shift, and the third light Z3 . Addi- 
tionally, if an O'Connel effect was visible in the light curves, we adjusted the spot 
parameters: its coordinates ((j),X), a spot radius r, and the ratio between the spot 
temperature and the surrounding surface. During the RV-curve modelling, we ad- 
justed the orbital semi-major axis a, the spectroscopic mass ratio q spec , and the 
mean system velocity v y . 

As the result of the iterative procedure of light- and RV-curve modelling, we ob- 
tained two groups of solutions. The first one comprises the evidently near-contact 
systems (V1010 Oph, WZ Cyg, VV Cet, DO Cas, FS Lup, and V747 Cen). The 
second comprises the systems whose light curves were better fitted by a contact 
configuration model, again with large temperature differences between the compo- 
nents (CX Vir, FT Lup, BV Eri, FO Hya, CN And, and BX And). 



8 



A. A. 




Fig. 1. Comparison between the observed and theoretical photometric (left panels) and RV (right 
panels) curves of V 1010 Oph, WZ Cyg, VV Cet and DO Cas. 
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Table5 

General information concerning spectral classification of our targets. 



system 


method 


spectral type 


author(s) 


WZCyg 


spectral class. 


F4V 


Ogloza & Rucinski (priv. comm.) 


DO Cas 


spectral class. 


A4V 


Ogloza & Rucinski (priv. comm.) 


BVEri 


spectral class. 


F2V 


Ogloza & Rucinski (priv. comm.) 


CNAnd 


spectral class. 


F5V 


Rucinski et al. (2000) 


BXAnd 


spectral class. 


F2V 


Belletal. (1990) 


CXVir 


spectral class. 


F5V 


Hilditch& King (1988) 


FTLup 


spectral class. 


F2V 


Hilditch et al. (1984) 


FO Hya 


B-V 


FOV 


Candy & Candy (1997), this work 


FS Lup 


B-V 


G2V 


this work 


V747 Cen 


B-V 


A5V 


Chambliss (1970) 


VV Cet 


B-V 


A5V 


Rahman (2000), this work 


VIOlOOph 


UV continuum 


A7V 


Corcoran et al. (1991) 





Fig. 2. Comparison between the observed and synthetic photometric (left panels) and RV (right 
panels) curves of FS Lup and V747 Cen. 
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Table6 

The results obtained from light- and RV-curve modelling of near-contact sytems. 
* - parameters fixed during light curve modelling (see Sections 3 and 4 for explanations). 



parameter 


V1010 Oph 


WZ Cyg 


VV Cet 


DO Cas 


FS Lup 


V747 Cen 


4 J 


86.14(37) 


83.497(32) 


81.89(13) 


89.94(30) 


82.64(23) 


87.16(18) 


T prim r „, 


/ JUU 


DjjU 


51 JU 


ojjU 


JoOU 


61 JU 


T S B C \K] 


5132(11) 


4932(2) 


6252(7) 


4297(26) 


5130(5) 


4275(16) 


^prim 


2.81773* 


3.1200(1) 


2.43060* 


2.4879(15) 


2.8181(5) 


2.51001* 


^sec 


2.8352(33) 


3.1200(1) 


2.4802(12) 


2.4751(1) 


2.8443(39) 


2.5342(14) 


<7spec 


0.470(3)* 


0.631(36)* 


0.284(8)* 


0.304(5)* 


0.470(7)* 


0.320(5) 


Vy[km/s] 


-19.65(4) 


1.05(88) 


-1.36(90) 


26.22(30) 


11.64(63) 


-18.55(46) 


a[R@] 


4.487(15) 


4.039(27) 


3.006(41) 


4.587(32) 


2.746(16) 


3.522(21) 


Lib 


11.323(12) 


10.4347(46) 


11.5776(67) 


12.3042(92) 


10.227(15) 


12.5055(64) 


L\v 


10.960(14) 


10.1042(44) 


11.2295(70) 


12.1172(99) 


9.978(16) 


12.2276(60) 






9.9063(42) 






9.805(16) 


11.9661(64) 


Ui 




9.4824(41) 






9.586(16) 


11.7091(67) 


Lib 


0.849(11) 


1.5717(38) 


1.0013(59) 


0.1278(64) 


2.404(16) 


0.1346(43) 


Liv 


1.122(13) 


1.8488(39) 


1.2078(63) 


0.2092(88) 


2.589(15) 


0.2156(58) 






2.1257(40) 






2.717(15) 


0.3034(72) 


Li, 




2.5910(38) 






2.918(15) 


0.4866(91) 




90.0* 




105.20(87) 




93.7(2.0) 


90.0* 




79.5(1.5) 




0.4(1.5) 




1.80(12) 


41.0(1.2) 


r[°] 


19.7(1.4) 




14.744(97) 




58.07(54) 


64.7(2.0) 


J spot 


1.378(31) 










1.2094(40) 


rj, irpprim 
1 spot /l e ff 






0.253(79) 




0.9136(10) 





4. Near-contact binary stars 

4.1. V1010 Oph 

For light-curve modelling, we used the only available BV data obtained by Le- 
ung (1974). As an O'Connel effect is visible in the light curves, we also introduced 
a spot into the model. During computations, we observed good convergence of 
Q.i to a critical value of the first Lagrangian point L\ , and finally a semi-detached 
configuration was obtained as the best solution. This result naturally explains both 
the presence of the hot spot on the secondary star that is present in our model and 
the period shortening seen in the O-C diagram (Kreiner et al., 2001), as the result 
of mass transfer from the more massive primary to the less massive secondary Al- 
though BFs (Fig. 3) were calculated from high signal-to-noise spectra, they did 
not reveal any features that could be interpreted as a hot spot. This may be a 
consequence of exceptional characteristics of this hypothetical spot in RV space, 
preventing it from being easily detected; or of insufficient spectral resolution; or of 
both these factors. 
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RV [km/s] RV [km/s] 

Fig. 3. BFs of V1010 Oph obtained with F2V RV standard in both quadratures. The phases calculated 
for the times of mid-exposure are given in the upper-right corner. 



4.2. WZCyg 

We took a new set of light and RV curves over a two-month period. Despite 
many attempts, we were not able either to remove or to explain the significant 
differences among the values of v Y determined separately for the RV curves of each 
component. High-resolution, good-SN spectra might resolve the ambiguity. Upon 
assuming a spectroscopic mass ratio of ^ spec = 0.631(36), we obtained a near- 
contact configuration (Table 6, Fig. 1) for this system. Because of the rather large 
uncertainty in q spec , we made two additional light-curve models with two extreme 
possible values of ^ spec (0.595 and 0.667), assumed constant in the computations. 
These models, however, likewise yielded a near-contact configuration. 

4.3. WCet 

As the data obtained with the R filter turned out to be affected by instrumental 
problems, we used only the BV data during light curve modelling. As a result, 
we obtained a semi-detached configuration, with the primary component filling its 
Roche lobe. Only by placing a cold spot on the surface of the primary component 
was it possible to obtain a good fit to the observed light curves. Since the primary 
component has a radiative atmosphere, this result indicates the possibility that the 
primary is of type Ap. Unfortunately, our BF profiles are not of a quality high 
enough for us to detect and confirm the existence of such a spot on the star. 

4.4. DO Cas 

As our three attempts to obtain new high-quality light curves failed, we used 
already-published BV data (Gleim & Winkler, 1969). Our light curve modelling 
used the spectroscopic mass ratio determined for the first time (<7 spec = 0.304) and 
resulted in a near-contact configuration. 
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4.5. FSLup 



We decided to improve the spectral classfication of FS Lup using our new data: 
a comparison of the spectra obtained at SAAO in the 3920-4750A region with 
respective spectra of RV standards indicated that the spectral type of FS Lup should 
be an early G. Additionally, from the B — V colour index, determined from our new 
photometric observations at the secondary minimum, we estimated the spectral type 
of the primary component as G2V (Table 5). The uncertainty deriving from the 
interstellar-reddening correction should not affect the effective temperature of the 
primary component (5860 K) by more than 150 K, as in the case of other systems 
with determined spectral types. 

To account for the O'Connel effect, a cold spot on the primary component was 
placed in that model which best reproduced the observed light curves. BF profiles 
obtained from medium-resolution spectra seem to reveal a cold spot on the primary- 
star hemisphere facing the secondary star, at coordinates in accordance with those 
obtained from light-curve modelling. However, spectra of higher resolution are 
necessary to confirm this finding. In investigating configuration, we found FS Lup 
to be a near-contact system, possibly semi-detached, with the primary star almost 
or entirely filling its Roche lobe (Table 6). 



4.6. V747 Cen 



The most interesting feature of V747 Cen is one which is, because of the high 
quality of our BFs, for the first time noticed in close binary systems: there is a well- 
defined region of increased intensity, visible in almost all profiles of the secondary 
component during second quadrature (Figs. 4, 5). We note that BFs of the sec- 
ondary component in first quadrature are also higher in the vicinity of the primary 
component (as compared with the corresponding BF behaviour of V1010 Oph in 
Fig. 3). This brighter region is underlined only in the plot obtained at phase 0.190, 
but one can easily locate it on the remaining plots presented in the left column 
of Fig. 4. The region may be a result of a hypothetical stream of matter, leaving 
the primary star through the first Lagrangian point, undergoing a slight deflection 
through centrifugal force, and then directly striking the side of the secondary-star 
photosphere visible in first quadrature. 

Light- and RV-curve modelling supported these conclusions: upon assuming 
that a hot spot on the secondary component is responsible for the O'Connel effect, 
we obtained a semi-detached configuration, with the primary component filling its 
Roche lobe (Table 6, Fig. 2). 
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Fig. 4. BFs of V747 Cen obtained with G5V RV standard. The brighter region on the secondary 
component in the vicinity of the primary component is underlined at phase 0.190 (see Sec. 4.6). The 
well defined brighter region visible in second quadrature (right panels) is indicated by an arrow in the 
first plot. Only one BF profile of the secondary component, obtained at phase 0.837, seems to be free 
from accretion effects. 
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Fig. 5. BFs of V747 Cen (as in Fig. 4) in phase-velocity space, rebinned in phase with a constant 
increment of 0.03. The bright region on the secondary component visible in second quadrature is 
marked. As expected, the profile of the radiative A5V primary component does not reveal dark spots. 



5. Contact systems with large AT 



5.1. CXVir 

In the high-quality BFs presented in Fig. 8, we notice even stronger differ- 
ences between secondary-component velocity profiles obtained in first and second 
quadrature than for the semi-detached system V747 Cen (Sec. 4.6). At phases be- 
tween 0.6 and 0.7, the triangular bright region (indicated by an arrow on the BF 
obtained at phase 0.747 in Fig. 8) dominates the signal of the secondary compo- 
nent. The latter begins to be distinguished from the brighter area in the BF obtained 
in phase 0.747 (Fig. 8, doubly underlined area on the same plot). The bright re- 
gion becomes narrower and fainter with phase, but remains visible even at phase 
0.844. BFs of the secondary component obtained in first quadrature are more in- 
tense (brighter) in the vicinity of the primary component (underlined area at phase 
0.202), exactly as with V747 Cen. This similarly suggests an ongoing mass trans- 
fer from the primary to the secondary. In Fig. 9, a few cold spots (i.e., the darker 
coherent trails) on the primrary are visible as well, outside the well defined bright 
region on the secondary star (indicated by longer marks). 

As the bright region on the secondary component introduces the largest effects 
in the BFs, we assumed one circular hot spot on the equator of the secondary star 
during light-curve modelling. We expected to obtain a semi-detached configuration 
as in the case of V747 Cen, but in fact we obtained a contact configuration, with a 
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Fig. 6. Comparison between the observed and synthetic photometric (left panels) and RV (right 
panels) curves of CX Vir, FT Lup, BV Eri and FO Hya. 
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Table7 

The results obtained from light- and RV-curve modelling of contact systems with large temperature 
differencies between components (CLdT). 
* - parameters assumed during light curve modelling (see Sections 3 and 5 for explanation) 



parameter 


CX Vir 


FTLup 


BVEri 


FOHya 


CNAnd 


BX And 




83.33(17) 


89.26(51) 


80.31(11) 


79.00(34) 


69.416(92) 


75.862(33) 


T e p f T[K] 

e JJ 


6450* 


6700* 


6700* 


7000* 


6450* 


6650* 




4694(26) 


4651(4) 


5387(20) 


4667(27) 


4726(34) 


4758(5) 




2 5146(8) 


2 723 in) 


2 3737(7) 


2 2714(25) 


2 6403(12) 


2 7760(5) 


Qspec 


343(3)* 


440(3)* 


0.274(2)* 


238(6)* 


385(5)* 


455(8)* 


vy\ Pint n / j 


-1 1 65(32) 


-2 30(30) 


-39 62(35) 


S4 42(30) 


-23 47(63) 


-24 00(60) 


a[Ro] 


4.584(19) 


3.537(13) 


3.365(12) 


3.184(21) 


3.163(22) 


4.424(36) 


f\%] 


21.2 


13.4 


20.0 


36.0 


3.0 


4.5 


L\u 






8.798(40) 








L\B 


9.350(51) 


10.597(23) 


8.920(40) 


10.746(81) 


9.952(18) 


11.0925(51) 


L\v 


9.242(48) 


10.315(24) 


8.748(37) 


10.886(78) 


9.830(16) 


10.8293(52) 


L\R 


9.039(49) 


10.164(24) 




11.091(74) 


9.710(15) 


10.6043(56) 


Lu 


8.956(47) 


9.909(22) 




11.093(96) 


9.531(13) 


10.2546(58) 


Liu 






0.996(23) 








Lib 


0.648(26) 


0.7151(46) 


0.946(22) 


0.361(13) 


0.781(36) 


0.9154(52) 


Liv 


0.789(30) 


0.9002(56) 


1.093(23) 


0.497(16) 


0.946(40) 


1.1300(58) 


Lir 


0.923(32) 


1.0850(61) 




0.630(18) 


1.114(43) 


1.3325(61) 


Lu 


1.185(35) 


1.3902(69) 




0.759(19) 


1.395(44) 


1.6874(60) 


hu 






0.1407(30) 








Hb 


0.1060(35) 


0.0481(16) 


0.1409(30) 


0.0894(63) 






hv 


0.1048(34) 


0.0463(18) 


0.1398(29) 


0.0672(60) 






l 3R 


0.1073(37) 


0.0456(18) 




0.0374(58) 






hi 


0.0957(37) 


0.0386(16) 




0.0304(74) 






4»[°] 


90.0* 


90.0* 


90.0* 


90.0* 


90.0* 


90.0* 


*H 


7.61(29) 


19.66(43) 


11.83(35) 


64.9(3.8) 


11.49(36) 


3.00(0.45) 


r[°] 


104.1(2.1) 


65.3(1.3) 


107.9(1.8) 


102.9(4.4) 


113.9(1.9) 


19.53(0.91) 


T /'rsec 
ispot/ 1 e ff 


1.2051(59) 


1.2053(33) 


1.141(4) 


1.1015(55) 


1.2138(72) 


1.177(15) 



large temperature difference between components and a significant value of third 
light h , comprising about 10% of the total flux (Table 7). Hilditch & King (1988) 
did not account for a third light in their analysis, as it is not obvious from their data 
whether the secondary minimum is total or partial. This apparently small detail 
always, however, has a large impact on the final light-curve solution. 

5.2. FTLup 

As in the case of V747 Cen and CX Vir, the BF profiles of the secondary com- 
ponent obtained in second quadrature reveal a bright region. In first quadrature, the 
BF profiles of the secondary are brighter in the vicinity of the primary (Fig. 10; 
for detailed explanations, check particularly the plots obtained at phases 0.233 and 
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Fig. 7. Comparison between the observed and synthetic photometric and RV curves of CN And and 
BX And. 



0.710). In Fig. 11, we display the BFs in the phase-RV plane. Here one can clearly 
see not only the track of the bright region (hot spot) in second quadrature but also 
several tracks of cold spots on the primary in both quadratures. 

Accretion processes between components considerably disturb the RV curve of 
the secondary, with v y as obtained from the RV curve of the secondary 10 km/s 
smaller than the v y obtained from the apparently undisturbed RV curve of the pri- 
mary. We consequently determined the spectroscopic mass ratio of this system us- 
ing only the secondary-component RV measurements obtained at second quadra- 
ture (Fig. 6). At these phases (Fig. 10, right panels), the influence of the well 
defined bright region can be easily taken into account during the RV measurment 
process. 

With a hot spot on the secondary component, we expected to obtain a semi- 
detached configuration from light curve modelling. However, as in the case of 
CX Vir, our best solution converged at the contact configuration with large AT 
between components and a non-negligible (4%-to-5%) contribution from a third 
light (Table 7). 
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Fig. 8. BFs of CX Vir obtained with F5V RV standard. The brighter region on the secondary compo- 
nent in the vicinity of the primary component is underlined at phase 0.202 (see Section 5.1). The well 
defined very bright region visible at second quadrature (right panels) is indicated by an arrow on the 
BF obtained at phase 0.747, at which point it begins to be distinguishable from the apparently pure 
profile of the secondary component (doubly underlined on the same plot). The outstanding accuracy 
of the BF profiles is noteworthy here. 
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Fig. 9. BFs of CX Vir in phase-velocity representation, rebinned in phase with a constant increment 

of 0.0235. In addition to the very bright region on the secondary component seen in second 
quadrature (as indicated by longer marks), cold spots on the primary component are visible. (One of 
these, near first quadrature, is indicated by shorter marks). 

5.3. BVEri 

The BF profiles of BV Eri again reveal a well defined brigter region on the sec- 
ondary component (Fig. 12), although not as clearly as in three systems previously 
discussed. Similarly, in first quadrature the BF profiles of the secondary are also 
brighter in the vicinity of the primary-component profile. 

To perform light-curve modelling, we used the only available very good-quality 
UBV data, published by Badee et al. (1983). In the model, we assumed a hot spot 
on surface of the secondary component responsible for the O'Connel effect and 
for the features discovered in the BFs. Our computations resulted in a contact 
configuration with large AT and with a significant (13% to 14%) contribution from 
a third light. 

5.4. FO Hya 

Our new B V (RI) C data confirmed the F0V spectral type of the primary com- 
ponent determined by Candy & Candy (1997). First RV curves yielded a con- 
siderably smaller value of the mass ratio (^ spe c = 0.238, in contrast with the value 
<7phot = 0.552 which Candy & Candy determined). Despite the relatively poor qual- 
ity of the spectra obtained at LCO, we additionally notice differences between the 
BF profiles of the secondary component obtained in both quadratures: the profiles 
appear brighter and are well defined in first quadrature, being considerably weaker 



20 



A. A. 




Fig. 10. BFs of FT Lup obtained by means of F2V RV standard. The region of increased intensity 
on the secondary component in the vicinity of the primary component is underlined in phase 0.233 
(see Section 5.2). The well defined high-intensity region visible in second quadrature (right panels) 
is indicated by an arrow on the BF obtained at phase 0.710, with the apparently pure profile of the 
secondary component doubly underlined on the same plot. 
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Fig. 11. BFs of FT Lup in phase- velocity representation, rebinned in phase with a constant 
increment of 0.025. In addition to the well defined region of increased intensity on the secondary 
component in second quadrature (indicated by longer marks), a few cold spots on the primary are 
visible. (One of these is indicated by shorter marks in first quadrature.) 




0.0 0.2 0.4 0.6 0.8 

phase 



Fig. 12. BFs of BV Eri in phase-velocity representation, rebinned in phase with a constant 
increment of 0.02. The well defined bright region of increased intensity on the secondary 
component (indicated by marks) is visible at phases in second quadrature. 
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in second quadrature. The profiles obtained before phase 0.75 seem to reveal a 
region of increased intensity, as with the four systems discussed above. 

Because the incompletness of the new light curve deprived us of evidence for 
an O'Connel effect, we first attemped to model the light curves without spots. This 
first model, however, did not yield a good fit. In our second attempt, using in- 
formation obtained from BF profiles and following Candy & Candy (1997), we 
introduced a hot spot on the secondary. The model from these computations, while 
yielding a reasonably good fit (Fig. 6), additionally delivered a contact configu- 
ration, a large AT between the components, and a consirable contribution from a 
third light (Table 7). 

5.5. CNAnd 

CN And is a very frequently observed and analyzed system, and different con- 
figurations were obtained from its light and RV curve modelling (Table 1). As 
almost all available archival light curves of CN And are of low quality, we decided 
to obtain new high-precision, multicolour (BVRI) data. Radial velocity curves were 
already obtained at DDO, and the spectroscopic mass ratio was derived by the BF 
method (Rucinski et al. 2000). In Fig. 13 we display binned BFs of CN And in the 
phase-velocity plane. The BFs, although individually rather noisy, give informa- 
tion of better quality when stacked, revealing a subtle track of a brighter region on 
the secondary component at second quadrature. We therefore decided to explain 
the O'Connel effect apparent in light curves by assuming one large hot spot on the 
secondary component. We obtained a marginal-contact configuration with large 
AT . However, in contrast with CX Vir, FT Lup, B V Eri, and FO Hya, it was not 
necessary to include a third light to obtain a good fit. 

We note that Zola et al. (2005) obtained a semi-detached configuration. Their 
model, however, assumed two cold spots on the primary component, whereas in 
our model a hot spot was placed on the surface of the secondary, in agreement with 
information obtained from BFs. 

5.6. BXAnd 

At first glance, our new BVRI light curves appeared to be free from an O'Connel 
effect. However, synthetic light curves obtained within a spotless model (yielding 
a contact configuration with / = 10% and large AT between components) did not 
fit the observed curves very well, with the discrepancy becoming especially pro- 
nounced around the secondary minimum. A closer inspection of the region around 
phase 0.3 reveals that it is slightly brighter than the one around phase 0.7. The 
assumption of a hot spot on the equator of the secondary component yielded a 
marginal-contact configuration with / = 4.5% and large AT (Table 7), but with 
about 40% smaller %;r ed . 

In order to validate the assumption of a hot spot in our model, we examined the 
available BF profiles obtained from medium-resolution spectra, but they did not 
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show any visible tracks of a hot spot on the secondary component. This outcome 
may be due to the small size of the spot (as obtained from light curve modelling), 
with detection perhaps only possible through BFs derived from high-resolution 
spectra. 

Table8 

Physical parameters of components derived from light- and RV-curve modelling. 



system 


Mi [M G ] 


M 2 [Mq] 


*i [*<=>] 


*2[*©] 


Tl[K] 


Ti[K] 


Li[Lq] 


L 2 [Lq] 


V1010 Oph 


1.887(19) 


0.887(15) 


2.01(3) 


1.40 


7500 


5132 


11.45 


1.22 


WZCyg 


1.589(32) 


1.003(37) 


1.70(4) 


1.37 


6530 


4932 


4.71 


1.00 


VV Cet 


1.042(43) 


0.296(20) 


1.48(5) 


0.79 


8150 


6252 


8.66 


0.85 


DO Cas 


2.130(45) 


0.650(25) 


2.22(4) 


1.29 


8350 


4297 


21.47 


0.51 


FS Lup 


1.301(23) 


0.611(20) 


1.23(3) 


0.85 


5860 


5130 


1.60 


0.45 


V747 Cen 


1.541(28) 


0.493(17) 


1.70(3) 


0.98 


8150 


4275 


11.42 


0.29 


CX Vir 


1.731(21) 


0.594(13) 


2.24(3) 


1.39 


6450 


4694 


7.78 


0.84 


FT Lup 


1.871(21) 


0.823(15) 


1.64(2) 


1.13 


6700 


4651 


4.86 


0.54 


BVEri 


1.560(17) 


0.430(07) 


1.70(2) 


0.96 


6700 


5387 


5.22 


0.70 


FOHya 


1.588(61) 


0.378(34) 


1.67(3) 


0.89 


7000 


4667 


5.68 


0.33 


CN And 


1.433(30) 


0.552(20) 


1.48(3) 


0.95 


6450 


4726 


3.40 


0.40 


BX And 


2.148(52) 


0.977(41) 


2.01(5) 


1.40 


6650 


4758 


7.08 


0.90 
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6. Physical parameters of components 



Using the results obtained from light- and RV-curve modelling (Table 6 and 
Table 7), we computed physical parameters of components (Table 8). Because 
the highly distorted components of contact binary stars significantly differ from a 
sphere, their effective radii R were calculated under the assumption that the surface 
area of each component is equal to the surface area of a spherical star of radius R . 
As the radii are scaled by their respective distances a between centres of masses 
of components, a -errors propagate accordingly on the radii of both components. 
Effective temperatures of primary components derived from spectral classification 
are known with an accuracy no better than 100-200 K, and these external errors, 
much larger than those obtained from light-curve modelling alone (Table 6 and 
Table 7), affect also the effective temperatures of secondary components. These 
uncertainties make it even more difficult to estimate errors in the luminosity L, but 
L errors are usually equal to 10% to 30% of L values (Yakut & Eggleton, 2005). 

Zola et al. (2006) considered a sample of contact binary systems for which 
physical parameters were derived from combined light- and RV-curve modelling. 
From the available data, they constructed mass-radius (M-R), mass-luminosity 
luminosity (M-L) and temperature-luminosity (T-L) diagrams. The M-L and T- 
L diagrams were also constructed using component physical parameters corrected 
for the effect caused by energy transfer from the primary to the secondary, by ap- 
plying equations from Mochnacki (1981), with an ad-hoc modification taking into 
account temperature differences between components (Zola et al., 2006). Theoreti- 
cal parameters of stars on ZAMS, TAMS, and hydrogen depletion in the core were 
computed with the EZ program (Paxton, 2005), which is based on the most recent 
version of Eggleton's STARS code (Pols et al., 1995). 

We have appended to the above-mentioned diagrams the physical parameters 
of the systems analyzed in this work (Table 8). As can be seen in the M-R dia- 
gram (Fig. 14), almost all primary components of our targets (the only exception 
is VV Cet) have parameters corresponding to the MS stars, while secondary com- 
ponents are oversized by a factor similar to those of components of W UMa-type 
systems. 

In Figure 15, we present two T-L diagrams. The first is constructed with the 
physical parameters from Table 8, while the second uses parameters corrected for 
an effect caused by an energy transfer between components. While the uncorrected 
parameters of the CLdTs secondary components (red squares) do not differ from 
those of the secondary components of near-contact systems taken from literature, 
in the second plot they show an obvious tendency to occupy the region avoided 
by secondary components of W UMa type systems. This matter will be discussed 
more thoroughly in the next section. 
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Fig. 14. The M-R diagram for systems analyzed in this work on the background of the sample of 
W UMa-type stars from Zola et al. (2006). The continuous line represents the theoretical lines of 
ZAMS and TAMS computed for metallicities of z = 0.01 and z = 0.02. "Near-contact" denotes the 
systems described in Sections 4.1-4.6, and "CLdT" the systems described in Sections 5.1-5.6. 
Physical parameters of these systems are listed in Table 8. 
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Fig. 15. T-L diagrams for systems analyzed in this work, without (left panel) and with (right panel) 
provision for effects of energy transfer. Parameters of the W UMa-type stars in the sample from 
Zola et al. (2006) are also shown. The second plot indicates that secondary components of CLdT 
systems lie in the region which secondary components of W UMa-type systems tend to avoid. 
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7. Discussion and conclusions 

We performed light- and RV-curve modelling of twelve close binary stars, for 
which contact configurations with large temperature differences between compo- 
nents have been reported at least once in the literature. In our analysis we applied 
new spectroscopic mass ratios, obtained through the BF method, determined for 
the first time for all systems but CN And. 

We obtained a near-contact configuration for V1010 Oph, WZ Cyg, VV Cet, 
DO Cas, FS Lup, and V747 Cen. The primary components of V1010 Oph, WZ Cyg, 
VV Cet, FS Lup, and V747 Cen fill (or nearly fill) their Roche lobes, and accord- 
ingly their secondaries are considerably oversized; an ongoing or past mass trans- 
fer can explain their larger radii according to the mechanism proposed by Webbink 
(1976) and Sarna & Fedorova (1989). In this configuration, an ongoing mass trans- 
fer should cause period shortening, but this is only confirmed in the O-C diagram 
of V1010 Oph (Kreiner et al., 2001). 

Among systems from this group, the most interesting results were obtained for 
V747 Cen and FS Lup. The latter appears to be the second-known short-period 
period (0.38 d) semi-detached system, which next to V361 Lyr (Hilditch et al, 
1997), is captured in a broken-contact stage of oscillation, as predicted by the TRO 
theory (see Section 1), or during the phase of a first-time approach to a contact 
phase. In the high-quality BFs (i.e., stellar images in velocity space) of V747 Cen, 
we discovered bright regions on the secondary component. The discovery is a sign 
of one or more hot spots, most probably a result of accretion occuring between 
the stars. This picture is in full accordance with the result obtained from light- 
and RV-curve modelling, which resulted in a semi-detached configuration for this 
system. 

From the light- and RV-curve modelling of CX Vir, FT Lup, BV Eri, FO Hya, 
CN And, and BX And, we obtained contact models with large temperature differ- 
ences between components (AT > 1000 K), in contradiction to theoretical predic- 
tions for contact systems. However, as in the case of V747 Cen, the most interesting 
results for CX Vir, FT Lup, BV Eri, FO Hya, and CN And were obtained from dis- 
covery of hot spots in BFs of their secondary components. As the existence of the 
hot spots can be explained only by accretion effects occuring between components, 
that discovery strongly supports semi-detached rather than contact configurations 
for these systems. 

The just-stated direct result aside, there additionally are indirect indications of 
the possibility that (at least some) contact models with large A T , obtained from 
light- and RV-curve modelling, are wrong. 

(1) From light-curve modelling of four systems with the strongest bright re- 
gions, we also obtained significant values of third light (Sections 5.1-5.4, and Ta- 
ble 7). The hypothetical third component is neither visible in BFs (as a slowly 
rotating sharp-line star) nor as a nearby field star in DSS images of these systems. 
However, without adding a third light into the model as adjusted parameter, in all 
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cases flat secondary minima and depth of primary minima were not reproduced 
well by these models, and as a consequence the obtained %^ ed weigh were at least 
a few dozens of percent higher. It is also worth noting that the bright regions are 
common features not only for CX Vir, FT Lup, BV Eri, FO Hya and CN And, but 
also for V747 Cen, for which light- and RV-curve modelling gave a semi-detached 
configuration without any third light, in full consistency with all the information 
inferred from its BFs. On this empirical basis, one may deduce that all the sys- 
tems are in fact semi-detached, with the duality of configurations obtained within a 
pure Roche-lobe model resulting from accretion processes taking place on different 
levels in the various individual systems. 

(2) The BFs of secondary components obviously appear different in the two 
quadratures. Whereas in first quadrature they do possess a characteristic smoothly 
increasing intensity toward primary components (this can be explained as a result 
of star-stream interaction producing a hot spot and causing considerable velocity 
turbulence), the well defined bright regions in second quadratures can be produced 
by returning streams encircling secondary components. If this is a correct picture, 
then these systems must be semi-detached. The positive (in absolute sense) depar- 
tures of RVs of secondary components from the Roche model, observed in the case 
of V1010 Oph, DO Cas (Fig. 1, right panel) and BX And (Fig. 7, right panel), in 
phases close to both minima, may be a sign of such high-velocity streams. Unfor- 
tunately, other systems were observed only during phases close to quadratures. 

(3) According to the O-C diagrams (Kreiner et al., 2001), continuous period 
shortening of FT Lup and CN And argues for a semi-detached configuration allow- 
ing for mass transfer from primary to secondary components. This empirical fact 
is consistent with information extracted from BFs of these systems, but it is incon- 
sistent with configurations obtained from light-curve modelling. CX Vir, BV Eri, 
FO Hya, and BX And do not reveal any obvious signs of a permanent period de- 
crease or increase. The latter is expected to be easily observable in O-C diagrams 
of contact systems with large temperature differences, as the mass transfer rate 
from secondary to primary star in common envelope is predicted to be of order 
lO~ 8 M /yr (Flannery (1976), Robertson & Eggleton (1977)). The lack of per- 
manent period increase in this group of systems is again in contradiction to the 
supposition of contact configurations. 

(4) The position on the T-L diagram of secondary components of contact sys- 
tems with large A T , constructed from their physical parameters as corrected for 
effects of energy transfer in a common envelope (Fig. 15, right panel) indicates that 
these stars differ significantly from the components of typical W UMa-type sys- 
tems. If their input physical parameters (Table 8) obtained within a false contact 
model from the WD code are wrong (due to the problems with unexplained sources 
of Z3), the corrections for energy transfer may in turn lead to very different output 
values, as observed in Figure 15 (right panel). 

(5) The discovery of bright regions in BFs of secondary components, the ob- 
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servation of significant differences in BF profiles in the two quadratures, and the 
discovery of spotted primaries in (at least) CX Vir and FT Lup (as described in Sec- 
tions 5.1-5.5) suggests that it may not be possible to derive accurate configurations 
from light-curve modelling. This is because smoothly varying EW- and EB-type 
light curves contain only information about changes of integrated binary-star light 
with phase, offering constraints too weak to yield a model consistent with the whole 
range of information extracted from spectroscopy. The modelling software based 
on the Roche model cannot properly handle all these effects, because there remains 
a large set of free parameters: our attempts to obtain much more reliable models of 
CX Vir, FT Lup, and CN And did not deliver unique fits. 

Although all the above findings strongly support the hypothesis of Rucinski 
(1986) and Eggleton (1996), that systems for which Roche-lobe-based computer 
models give contact configurations with A T > 1000 K are in fact semi-detached 
(see Section 1), further observations are necessary for full confirmation of this find- 
ing. Doppler-imaging observations, as made by Pribulla & Rucinski (2008) for 
AW UMa, may be needed. The cited observations have revealed that the secondary 
component has the aspect of a small core-like star surrounded by optically thick 
matter. Further, the recent discovery of a double-peaked H a line in a short-period 
(1.6 d) SV Cen (Siwak et al., 2009) and (0.42 d) TYC 2675-663-1 W UMa type 
systems (Caballero-Garcia et al., 2010) may indicate an existence of accretion disk 
around one component, and a semi-detached rather than a contact configuration. 
This is a possibility previously envisaged by Zola (1995). The echelle spectra of 
CX Vir and FT Lup, covering a wide spectral range (UV-IR), do not reveal emis- 
sion components in hydrogen lines. Further progress on contact-binary systems 
with large A T between components may consequently be obtained from mod- 
elling of the available and the prospective new BFs, obtained over a whole orbital 
period. 
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